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INTRODUCTION 


Reflectivity  is  associated  with  an  interface  between  media  of  different 
refractive  indices.  The  front  of  a  mass  density  shock  wave  causes  a  change  in 
refractive  index  by  forcing  a  change  in  the  spatial  density  or  character  of  the 
optical  scattering  sites  and/or  the  magnitude  of  electric  moment-field  interac¬ 
tion  matrix  elements.  It  has  thus  been  clear  for  some  time  that  employing  a 
propagating  shock  front  as  a  reflectivity  surface  permits  the  study  of  details  of 
the  transition  between  the  shocked  and  the  unshocked  states,  as  well  as  the  pos¬ 
sibility  of  understanding  the  molecular  level  contributions  to  the  index  of 
refraction. 

Until  now  shock  front  reflectivity  experiments  have  been  used  primarily  as  a 
tool  for  estimating  the  thickness  of  the  shock  transition  region  (refs  1  to  3), 
with  the  classic  work  being  that  of  Cowan  and  Hornig  (shock  front  thickness  in 
gases)  (ref  I).  More  recently,  however,  we  have  attempted  a  very  detailed  com¬ 
parison  between  theoretical  and  experimental  reflectivities  for  shocked  liquid 
water  at  5.8  kbar  (0.58  GPa)  and  have  found  a  difference  between  theory  and 
experiment  which  cannot  be  explained  by  employing  a  refined  (secondary  reflec¬ 
tion)  Fresnel  reflectivity  theory  (ref  4)  (fig.  1). 

A  large  discrepancy  exists  between  the  theoretically  predicted  and  experi¬ 
mental  reflectivities  for  the  parallel-to-the-plane-of-incidence  optical  polari¬ 
zation,  as  illustrated  in  figure  1.  That  discrepancy  increases  as  the  optical 
angle  of  incidence  increases  or,  equivalently,  as  the  direction  of  the  E  field  of 
the  reflecting  beam  becomes  more  and  more  parallel  to  the  direction  of  shock 
propagation  (see  insert  fig.  1).  Because  the  case  of  the  perpendicular-to-the- 
plane-of-incidence  polarization  yields  quite  good  agreement  between  theory  and 
experiment,  the  parallel-to-the-plane-of-incidence  dependence  on  angle  strongly 
suggests  that  the  shock  front  must  be  strongly  influencing  those  terms  in  the 
optical  scattering  matrix  which  have  an  initial  and/or  final  electric  dipole 
moment  state  in  the  direction  of  shock  propagation. 

An  analogous  case  exists  where  a  low  power  laser  probe  beam  is  used  to  sam¬ 
ple  a  material  which  has  had  optical  phonon  states  pumped  with  a  preferred  direc¬ 
tionality  by  the  imposition  of  a  high  powered  laser  beam  (stimulated  Raman  scat¬ 
tering)  (ref  5).  In  that  case,  the  interaction  of  the  probe  beam  and  the  optical 
phonon  diminishes  the  probe  beam  intensity  and  creates  side  band  beams  (Stokes 
and  anti-Stokes)  which  are  spatially  separated  from  the  probe  beam. 

This  report  demonstrates  that  the  one-dimensional  strain  shock  front  region 
(the  transition  region  between  two  differing  mass  density  hydrostatic  states  of 
water)  could  serve  as  a  mechanical  pump  in  the  creation  of  excited  (i.e.,  other¬ 
wise  relatively  unpopulated)  electric  dipole  states  in  the  liquid  water.  The 
reflectivity  measuring  laser  beam  then  would  become  equivalent  to  the  above- 
mentioned  probe  beam,  and  that  reflectivity  beam  would  interact  with  the  excited 
states  of  the  shock  front  region  forming  side  bands.  This  would  clearly  have  an 
effect  upon  the  measured  reflectivity.  For  example,  the  finite  angular  width  of 
the  photo  diode  means  that  the  non-interacted  part  of  the  reflectivity  beam  and  a 
spatially  separated  side  band  are  not  simultaneously  imaged  on  the  photo  diode. 
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SHOCK  FRONT  PUMPED  EXCITED  STATES 


The  probe  beam,  excited  state,  and  side  band  geometry  shown  in  figure  2 
assumes  that  the  excited  state  wave  vector,  k  ,  is  parallel  to  the  direction  of 
shock  propagation.  For  probe  (reflectivity)  $eam  frequency  oj  and  excited  state 
energy  the  conservation  of  momentum  and  energy  leads  to 

k'  -  k  +  k  (1) 


k" 

£  " 


(2) 


If  the  z  (shock  propagation  direction)  component  of  equation  1  is  combined  with 
equation  2  and 

k^  *  k'  cos  (  9  -  \Ji)  (3) 

z 


then 


(i) 

—  -  (l  +  — )  COS  (  0  -  1(0  -  cos  0 


(A) 


so  that  knowledge  of  0,  ^  and  yields  a  value  for  k  • 

By  manipulating  equations  1,  2,  and  3,  it  can  be  shown  that 


0) 

sin2  0  =*  -  (l  +—  f  sin2  0  cos2  + 


o>  „ 

+  2  (l  +  — )  sin  0  cos  0  sin  \J>  cos  \p  + 


U) 

-  (l  H - -  j2  sin2  cos2  0  + 


0) 

+  2  (l  +  — -)  sin2  0  cos  ip  + 


-  2  (l  +  -jjp)  sin  0  cos  0  sin 


or,  to  first  order  in  tp, 


sin2  0  *  -  (l  +  -~f  sin2  0  +  2  \J>  (l  +  sin  0  cos  0  + 

w  2  <*> 

+  2  (l  +  — )  sin  0-2  <J»  (l  +  — -)  sin  0  cos  0 


(5a) 


(5b) 
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Equation  5b  yields  for  »j> 


Employing  the  same  approach  as  above,  it  is  possible  to  show  that, 

makes  an  angle  $  with  respect  to  the  direction  of  shock  propagation,  then, 
to  first  order  in  ij>, 


(6a) 

if  k 
rv 
again 


Because 

imately 
will  be 


CO  8 

L  c 


6-  cos 
cos  9  sin  9 


*] 


(6b) 


the  difference  between  tan  0  and  r 


2  2 
cos  <j>  -  cos  9 


*  is  small  for  0  approx- 


cos  9  sin  0 

75  degrees  and  <(>  approximately  10  degrees,  only  given  by  equation  6a 
considered  in  the  following  analyses  (i.e.,  <p  *  0  is  assumed). 


Equation  6a  and  figure  3  show  that  w  >  0  (corresponding  to  anti-Stokes 
radiation)  gives  >  0  and  represents  the  earliest  signal  to  be  imaged  on  the 

photo-diode.  oj^  <  0  (corresponding  to  Stokes  radiation)  gives  ^  <  0  and  repre¬ 
sents  the  last  signal  (as  the  shock  front  sweeps  down  the  sample)  to  be  imaged  on 
the  photo-diode. 


The  essential  idea  is  that  within  the  shock  front  region,  optically  active 
excited  states  separated  by  energy  ho)  are  caused  to  exist  and  transitions 
between  those  states  make  up  reflectivity  side  bands  when  the  momentum  conserving 
k  states  are  available.  Normally  such  side  band  signals  are  observed  spectro¬ 
scopically  by  studying  optical  intensity  as  a  function  of  frequency.  In  the  case 
of  the  measurements  of  the  shock  front  optical  reflectivity  discussed  here,  how¬ 
ever,  it  is  the  spatial  separation  associated  with  the  side  band  signals  that  is 
potentially  observed. 

A  typical  parallel-to-the-plane-of-incidence  (optical  polarization)  noise 
averaged  photo-diode  record  is  illustrated  in  figure  4.  If  (tc  -  tg)  is  taken  as 
the  separation  between  the  main  reflection  signal  and  the  anti-Stoxes  signal  at 
the  converging  lens  photo-diode  system,  then 


(tc  -  V 


\{>  R  gin  9 


U 


S 


(7) 


where  Ug  is  the  shock  velocity,  and  R  is  the  distance  from  the  scattering  volume 
to  the  converging  lens  photo  diode  system.  Equation  7  assumes  scattering  volume 
dimensions  are  small  compared  to  R  (which  was  satisfied  experimentally),  and  that 
the  beam  diameter  is  less  than  beam-side  band  separation  at  the  photo-diode.  The 
laser  beam  always  had  a  diameter  of  approximately  0.3  mm  at  the  scattering  volume 
(prior  to  scattering).  It  is  assumed  that  the  main  (scattered)  reflectivity  beam 
was  less  than  1  mm  in  diameter  at  the  photo-diode. 
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For  the  typical  experimental  values  of  R  *  23  cm,  Ug  ~  2.1  x  lO'*  cm/sec  at 
5.8  kbar  (ref  6),  (t^  -  tg)  ~  0.8  y  sec,  and  optical  angle  of  incidence  of  78 
degrees 


i|>Ag  (78  degrees)  ~  0.43  degree 


(8a) 


0.27  degree  corresponds  to  ~  1  mm  at  23  cm.  Thus,  from  equation  6a, 

,-3 


G) 

r~3.2x  10 


(8b) 


For  a)  =  2xf  and  f  =  (5145  A)  *  corresponding  to  the  wavelength  of  the  reflectiv¬ 
ity  laser  employed  in  these  experiments,  equation  8b  givest 


f  ~  61  cm 
v 


-1 


Equation  8c  substituted  into  equation  4  yields 

IT  I AS  ~  8*°  x  10’3 


(8c) 

(8d) 
-1 


If  the  value  for  k  within  the  scattering  medium  (water)  is  used,  then  k  =  2xnX 
with  n  =  1.33  (the  index  of  refraction)  and  X  =  5145  A  gives  from  equation  8d 


k  | 

\i  1 


AS 


(78  degrees)  ~  1.3  x  10^  cm  * 


-1 


<8e) 


The  f  ^  61  cm  value  from  equation  8c,  which  actually  is  good  to  only  one 
significant  figure,  is  to  be  compared  with  Raman  signals  observed  near  175  cm~J 
and  60  cm"*  by  Walrefen  in  water  (ref  7),  and  near  190  cm"~*  and  50  cm 
(inferred)  infrared  signals  observed  in  ice  by  Bertie  and  Whalley  (ref  8).  The  ~ 
near  60  cm“*  (50  cm~*)  signal  has  been  associated  with  transverse  displacements 
(bending)  of  the  intermolecular  OH  bond  in  tetrahedrally  bonded  water  (refs  7  and 
8). 


0) 


Setting  (— )  ^  -  3.2  x  10  ^  in  equation  8b  yields  the  corresponding  Stokes 
a  corresponding  separation  angle  (again  as  a  first  order  quantity) 


line  with 
given  by 


(78  degrees) 


*  0.43  degree 


(9a) 


t  Spectroscopic  units  (i.e.,  cm  )  are  used  for  frequencies  in  this  report. 
Ordinary  frequency  units  are  obtained  by  multiplying  the  spectroscopic  units  by 
the  speed  of  light. 
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with  the  corresponding  result 


k 

v 

IT 


-  8.0  x  10 


-3 


(9b) 


and 


V  Is  * 


I  .3 


,  3  -1 

10  cm 


(9c) 


The  minus  sign 
conserving  moineni 
opposite  to  the 

parallel  to  the 


<atlon  9c  icdtdt^  that  the  phun  Oil  *  »te  respon^lt  for 
*  the  Stokes  process  is  4n  emitted  pUoKiti  which  is  «i  1  tel 
; t  i on  of  shock  prap4f!at  Ion,  Reafcjftbetvd  that  onLy  k 

■ction  of  shock  propagation  are  core?  i  tiered  here 


i* 


Q  in 


equation  6b  ha  a  bb  ,  assumed). 


DISCUSS IQN 


The  number  o«n  states  coat r ibatlfig  to  Che  7k  r,<  m  n  tl-v  ,I<<  *ig- 

nais  can  be  esc  l  \  uom  the  known  refle  ilvlty  cue  i  i  i  c  (  5%  at  /8  legrec* 

angle  of  incident  i*)r  the  1  watt  CW  laser  beam  which  it-  scattered  into  the 
lens-photo-diode  for  ms tley  7  -  it  it  is  assumed  that  the 

Stokes  and  anti-St  !i  i  signals  have  intensities  appr  ximately  equal  t<  the  si  ipie 
reflectivity  (which  tppeara  to  he  the  experimental  situation  at  5-8  kbar  and  78 
degree  angle-of-incidance)  then,  since  approximately  1  erg  oi  the  laser  beam  is 
converted  into  Stokes  ur  anti-Stokes  radiation, 


N  r  nuiihei  of  side  bant  photons 
P 


N  —  3  x  10**  photons  (  10) 

P 


The  comparis  f  the  N  value  with  the  expected  room  temperature  thermal 
equilibrium  phono  occupation  for  w  (k  )  states  is  interesting.  Whether  the 
expression 


*  "ho) 

N(1)  (T)  =  {  exp  (  ~ 

v  2ir  1  ^  v  kT 


)  -  1  J  1  dk^ 


Ola) 


is  employed  for  one-dimensional  physics  or 
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(3) 

v 


(T) 


V 

(2tt)3 


e*p  (  IT  >  ■ 


k2  d  k 
v  v 


(lib) 
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for  three-dimensional  physics,  with  dk  ~k  ,  and  whether  V  -  1  crn  ^ ,  L  *  1  cm  or 
expected  (ref  9)  shock  front  thickness  and  laser  beam  cross  section  related 
values  (V  «  2  x  10~^  cnr  ,  L  ■  20  A)  the  result  is  N  (T)  «  N  by  at  least  six 
orders  of  magnitude.  ^ 

Such  large  differences  between  N  (T)  and  Np  have  already  been  observed  (ref 
10)  for  laser  pumped  stimulated  Raman  spectroscopy  (SRS)  experiments  and  implied 
in  theoretical  considerations  by  Coffey  and  Toton  (ref  11). 

The  implication  of  the  difference  between  N  (T)  and  Np  is  clear.  The  shock 
front  region  pumps  up  phonon  states  to  occupancies  far  greater  than  either  the 
thermal  equilibrium  state  before  the  shock  front  or  the  thermal  equilibrium  state 
after  the  shock  front.  The  shock  front  region  should  thus  be  viewed  as  physi¬ 
cally  different,  in  the  strongest  sense  possible,  from  the  overall  entity  known 
as  a  shock  wave.  This,  if  correct,  represents  a  first  for  mass-density  shock 
wave  physics:  the  availability  of  an  experimental  tool  capable  of  observing  the 
states  within  the  shock  front  region  in  detail. 

Any  study  which  claims  that  mass  density  shock  pumped  Raman  spectroscopic 
effects  are  making  a  contribution  to  shock  front  reflectivity  measurements  must 
attempt  to  calculate  the  intensities  of  the  Raman  lines.  It  is  relatively  easy 
to  show  that  any  intensity  calculation  based  on  equilibrium  effects,  such  as  that 

contained  in  the  Rayleigh  ratio  value  R^  0  ~  10  ^  cm  *  for  water  (ref  12)  will 

give  intensities  at  the  photo-detection  system  approximately  six  orders-of- 
magnitude  smaller  than  are  required  to  explain  observations.  That  is  not  sur¬ 
prising  in  light  of  the  previous  comments  concerning  the  difference  in  phonon 
occupation  within  the  shock  front  region  as  compared  to  the  equilibrium  state. 
Fortunately,  it  is  possible  to  carry  out  a  simple  calculation  for  the  Raman  scat¬ 
tered  field  strength  based  upon  a  simply  stated  nonequilibrium  contribution. 

The  scattered  field  strength,  E"  ,  at  the  detection  system  may  be  shown  (ref 
13)  to  satisfy  the  relation 

- 2  TT2  I  G  1 2  sin2  ||> 

in2-  2*  2  4 —  <12> 

R  z  \ 
o 

where 

I  G~  I2  "  F-oi  E„k  /  ua.  dv  /  5e«k  dv  <13) 


RQ  is  the  distance  between  the  photo-detection  system  and  the  scattering  volume, 
e  is  the  scalar  dielectric  constant,  is  as  employed  in  figure  2,  E  is  the 
incident  field  strength,  is  the  instantaneous  amplitude  of  the  fluctuation 

of  the  permittivity,  dV  is  a  differential  element  of  the  scattering  volume,  and 
the  bar  indicates  an  average  over  the  molecular  motion. 


Let  us  approximate  the  contents  of  equation  13  by 


G 


i 


V 


scatt 


(14) 
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and,  appropriate  to  a  one-dimensional  (strain)  shock  problem, 


«e  i  (6c)  6U  6  (15) 

where  the  delta  function  subscript  1  indicates  the  direction  of  shock  propaga¬ 
tion.  Clearly  equation  15,  when  applied  to  the  perpendicular-to-the-plane-of- 
incidence  optical  polarization  case,  yields  ■  0.  The  parallel-to-plane-of- 
incidence  case,  where  E  is  perpendicular  to  k  and  in  the  plane  formed  by  k  and 
k-*,  yields 


G.  ,  -  e.  (6e)  E  V  sin  9 

11  1  o  scatt 


(16) 


where  is  a  unit  vector  in  the  direction  of  shock  propagation  and  9  is  the 
angle  of  optical  incidence. 

If  i  is  the  laser  beam  diameter,  then  consideration  of  the  projection  of  the 
beam  diameter  on  the  shock  front  within  the  plane  of  incidence  gives 


L  l 


scatt  cos  9 


(17) 


where  is  the  shock  front  thickness.  Combining  equations  12,  16,  and  17  results 
in 


2  2  4  2 

iE',2  it  L  1  _  2  „  ,  2  , 

|— | - 4  -  —  [— J  tan  9  sin  * 

o 


a4r  2 

O 


(18) 


Substituting  from  equation  6a  then  gives 

.  2  2  4  ~  2  a)  2 

iE  i2  it  L  l  f6e 


E  ,,4  2  '•e 

o  4A  R 

o 

where  «  w  has  been  assumed. 


/■OEi  r  \)\  4  „ 

f — '  ( —  tan  9 

'•e  '  v'u>  ; 


(19) 
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Substituting  L  ■  20  A,  l  *  1  mm,  A  -  5  x  10  cm,  and  RQ  ■  23  cm  into  equa¬ 
tion  19  yields 


||1|2  -  3.0  x  10"3  (|£)  tan4  9 


(20a) 


If  (— )  ■  1,  representing  a  strong  non^equilibrium  shock  front  induced  permit¬ 
tivity  fluctuation,  is  used  along  with  (— )  ~  3.2  x  10_3  from  equation  8b  and  0  ■ 
78  degrees  the  result  is 
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(20b) 


|’!“|2  ~  1.5  x  10  5 


While  the  reflectivity  contribution  indicated  by  equation  20b  is  approximately 
three  orders-of-magnitude  smaller  than  necessary  to  explain  experimental  observa¬ 
tions,  it  is  at  the  same  time  approximately  three  orders-of-magnitude  larger  than 
a  typical  equilibrium  Rayleigh  ratio  contribution. 

The  sin  term,  with  small  in  equation  18  i3  the  principal  reason  for 
equation  20b  being  smaller  than  necessary.  There  is,  however,  another  Stokes  and 
anti-Stokes  contribution  which  does  not  suffer  from  the  small  if/  effect.  To  this 
point  we  have  only  considered  the  Stokes  and  anti-Stokes  contributions  relative 
to  (generated  by)  the  simple  reflected  beam  labeled  k  in  figure  2.  It  is  also 
possible  Co  have  a  contribution  from  the  same  scattering  volume  relative  to  (gen¬ 
erated  by)  the  incoming  laser  beam. 


For  the  one-dimensional  permittivity  fluctuation  of  equation  15  it  is  possi¬ 
ble  to  show  that,  for  an  incoming  beam  k  making  an  angle  0  with  respect  to  the 
normal  to  the  shock  front,  equation  18  becomes 


o 


2 


2  2 
L  i 


2  2 
tan  0  sin 


(21a) 


where  again  E'  is  non-zero  only  for  the  parallel-to-the-plane-of-incidence  opti¬ 
cal  polarization,  and  Z'  itself  is  entirely  within  the  plane  of  incidence. 
Employing  the  same  values  for  the  parameters  in  equation  21  as  were  employed  with 
equations  20  gives 


||n2  -  5.3  x  10  2  (21b) 


An  approximately  5%  Stokes  and  anti-Stokes  reflectivity  contribution  is  more 
than  sufficient  to  explain  the  larger  than  expected  water  reflectivity  signals 
for  the  parallel  to  optical  plane  of  incidence  case  (ref  4). 

The  effect  associated  with  equation  21a  also  has  the  desirable  properties  of 
an  increasing  contribution  for  an  increasing  angle  of  incidence  0,  and  a  zero 

contribution  for  E  perpendicular  to  optical  plane  of  incidence.  Both  properties 
^o 

are  in  agreement  with  experimental  observation  (ref  4). 

The  geometry  associated  with  equation  15,  which  is  equivalent  as  a  phenome¬ 
non  to  saying  that  k^  is  normal  to  the  shock  front,  restricts  E'  in  equations  18 

and  21a  to  be  entirely  within  the  plane  of  incidence.  If  the  equation  15  condi¬ 
tion  is  relaxed,  then  two  effects  occur: 
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1.  A  parallel-to-the-plane-of-incidence  incoming  beam  will  have  a  Stokes 
and  anti-Stokes  reflectivity  contribution  with  perpendicular-to-the-plane-of- 
incidence  components. 

2.  A  perpendicular-to-the-plane-of-incidence  incoming  beam  will  have 
non-zero  Stokes  and  anti-Stokes  reflectivity  contributions. 

The  first  of  these  two  effects  is  potentially  the  most  interesting,  because 
such  an  effect,  not  predicted  to  exist  for  simple  reflectivity,  represents  a 
decisive  experiment  which  can  be  simply  carried  out  in  a  typical  shock  wave  lab¬ 
oratory.  For  the  parallel-to-the-plane-of-incidence  reflectivity  experiment  a 
polarization  analyzer  at  the  photo-detection  system  must  be  added.  The  analyzer 
would  be  polarized  perpendicular  to  the  incidence  so  that  only  Stokes  and  anti- 
Stokes  contributions  would  be  observable. 

For  tetrahedral  hydrogen  bonded  water  there  will  be  some  shock  compressed 
hydrogen  bonds  at  non-zero  angle  relative  to  the  direction  of  shock  propagation. 
Thus  it  is  reasonable  to  expect  an  out  of  the  optical  plane  of  incidence  polari¬ 
zation  contribution,  although  perhaps  small  in  amplitude,  if  the  basic  ideas  con¬ 
tained  within  this  paper  are  correct. 


CONCLUSIONS 


We  have  postulated  Stokes  and  anti-Stokes  contributions  to  the  shock  front 
optical  reflectivity.  Those  contributions,  viewed  here  as  arising  only  from  the 
shock  front  region,  represent  a  new  and  important  tool  for  probing  the  transition 
between  the  unshocked  and  shocked  states. 

In  the  near  future  these  authors  intend  to  look  for  perpendicular-to-the 
plane-of-incidence  signals  with  a  parallel-to-the-plane-of-incidence  incoming 
beam.  Also,  employing  laser-generated  shock  waves  in  water  at  approximately 
5  kbar  amplitude,  the  Ultrafast  Spectroscopy  and  Laser  Physics  Laboratory  at  the 
City  College  of  New  York  will  soon  search  for  Raman  frequency  shifted  lines 
within  the  shock  front  reflected  signal. 
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Theory  for  5.8  kbar  from  reference  4.  Insert  shows  reflec¬ 
tivity  geometry  for  parallel-to-the-plane-of-incidence 
optical  polarization.  Reflectivity  beam  is  a  laser  at  5145 
A  (vacuum). 


Figure  1.  Water  reflectivity  data 
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/  IMAGING  AND 
/  PHOTO  DIODE 


k  is  the  main  reflectivity  beam  wave  vector,  k"  the  sideband 
wave  vector,  and  0  the  optical  angle  of  incidence. 


Figure  2.  Conservation  of  k  vector 


14 


dp  is  the  photo-diode  effective  diameter,  and  tj  and  t2 
represent,  respectively,  the  earliest  and  latest  times  for 
the  main  reflectivity  beam  to  be  imaged  on  the  photo-diode 
(for  vanishing  beam  diameter  d^) .  Us  is  the  shock  velocity. 


Figure  3.  Photo-diode  imaging  geometry 
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78  degree  angle  of  incidence  and  parallel-to-the-plane-of- 
incidence  optical  polarization. 


Figure  4.  Noise  averaged  photo-diode  record 
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